Radiative Corrections to the Off-Resonant Maxwell-Garnett Dielectric 
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We compute the dielectric constant of a Maxwell-Garnett (MG) dielectric renormalized by radiative correc- 
tions at off-resonance frequencies. To this aim, the self-energy of the dipole constituents of the dielectric 
is calculated following Ref. Q]. The spectrum of coherent emission is obtained in exact form for the MG 
model. Contrarily to the common assumption, its expression contains an only local field factor. Beyond the 
MG model, leading order corrections due to recurrent scattering are implemented. Several experiments are 
proposed to validate our results in both three and two dimensional samples. 



I. INTRODUCTION 

The effective dielectric constant in nonpolar fluids, 
spherical suspensions, metallic nanoparticle composites 
and generic disordered media is a problem profusely stud- 
ied along the history. Its study goes back in time to 
the seminal works of Claussius^, Mossotti^, Lorentai, 
Lorena^ and Maxwell-Garnett 6 . The problem has re- 
ceived the attention of a number of Physics commu- 
nities. These include Fluid and Statistical Physics^— , 
Chemical Physics^ - - — , Stellar Atmospheres and Inter- 
stellar Dus1?i£~— , Multiple-scattering Optics^r— , Quan- 
tum Optics and Fluorescence^ - —. The underlying fun- 
damental theory is that of linear transport in disordered 
media. Any analytical approach to the problem is based 
on some model and, even for the simplest ones, further 
approximations are needed to achieve meaningful expres- 
sions. In all the works cited above the modelling con- 
sists of treating the dielectric as a collection of identi- 
cal point dipoles which are distributed in space accord- 
ing to some equilibrium statistics. The construction of 
the dielectric constant is therefore a bottom-up approach 
in which knowledge about the relevant microscopical de- 
tails is assumed. The problem reduces this way to one 
of multiple-dipole interactions. Its solution involves the 
calculation of a stochastic kernel (in statistical terms) , an 
effective potential (in many-body terminology), an elec- 
trical susceptibility (in optics) or a self-energy function 
(in quantum terminology). To this aim, cluster expan- 
sions have been developedii&^i. They all base on per- 
turbative series. Thus, the problem turns into one of 
finding out the appropriate parameter of the expansion 
and of selecting the terms of the resultant series which are 
relevant in some plausible approximation. The authors 
of most of those works justify their analytical approxi- 
mations with numerical simulations. Numerics have the 
advantage that one can isolate the physical aspects of in- 
terest by just ignoring those aspects which are not being 
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studied. However, sometimes this leads to oversimplifica- 
tions which are not physically acceptable. For instance, 
the radiative emission of an excited dipole generally de- 
pends on its coupling to both electrostatic and radiative 
modes. We will argue that those approaches which ne- 
glect the contribution of longitudinal modes to radiative 
emission cannot yield satisfactory results. On the other 
hand, those approaches which incorporate only electro- 
static interactions cannot be consistent with the optical 
theorem. They give incomplete results for the line-width 
of the spectral function. 

A more technical aspect is that of the effect of the 
spatial dispersion of the dielectric function in the renor- 
malization of single-particle polarizabilities. The spatial 
dispersion and its effects depend on the specific micro- 
scopical structure of the medium. In general, further 
modelling is needed. Nevertheless, there are still aspects 
which can be considered generic under a minimum of re- 
strictions. It is the first aim of this article to study the 
necessary conditions for this to be the case. For this pur- 
pose, we analyze the role of the different length scales 
which enter the problem and identify the appropriate ex- 
pansion parameters. This allows us to compute, in good 
approximation, the effective dielectric constant, e e //, of 
a statistically homogeneous medium off-resonance. The 
computation at resonance is left for a separate paper as 
it involves a different expansion. 

Closely related to the computation of e e // is the cal- 
culation of the spectrum of stimulated emission. To this 
aim, we exploit our results in Ref. [lj. We address this 
issue in detail as it serves us to clarify some usual mis- 
conceptions regarding the role of local field factors in the 
spectrum of dipole emission. 

Our object of study is a complex dielectric medium 
made of point dipoles embedded in the free space of di- 
electric constant eo- The latter can be generalized to 
allow for the embedding in an auxiliary medium which 
behaves as a continuum background with well-defined 
passive constant permitivity. The spatial distribution of 
dipoles is statistically isotropic and homogeneous. The 
dipoles polarize the electromagnetic (EM) vacuum. The 
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polarized vacuum contains additional fluctuations with 
respect to (w.r.t.) the in- free-space vacuum. Those 
fluctuations are referred to as polaritons when coherent. 
As a result, both the coherent transport features of the 
medium and the self-energy of the dipole constituents 
get modified. In turn, that implies that both the dielec- 
tric constant and the single-particle polarizabilities, a, 
get renormalized. Our treatment rely on a dipolc-Born- 
Markov approximation as explained in Ref. [l| and is a 
Green's function based approach. 

In particular, we deal in this work with a Maxwell- 
Garnett dielectric. The reasons being that analytical 
results can be obtained, comparison with previous ap- 
proaches can be made and experimental tests can be 
proposed. Generically, we consider a Maxwell- Garnett 
dielectric as a random medium made of well-separated 
point dipoles of bare electrostatic polarizability a?o whose 
correlation length £ satisfies i( « 1 for the frequencies 
of interest, uj = ck. Well-separated means that the dis- 
tance between dipoles is a few times greater than their 
longitudinal dimension in order for the dipole approx- 
imation to be valid. In the long-wavelength limit of 
the effective theory, q£ — > 0, no additional assumption 
about the degree of correlation is needed and the rele- 
vant correlation function for long wavelength modes is 
the two-point function, /t(r) & 21 i 30 . It contains at least 
two inherent terms. Those are, one which stands for the 
exclusion volume around each dipole and another one 
which stands for the self-correlation. Further, if the self- 
correlation term is discarded, no matter the exact form 
of the exclusion volume nor the rest of the terms in h(r), 
the resultant model obeys the usual MG formula for the 
effective susceptibility^ ' 23 ' 31 (see below). In this work 
we will refer to this approximation as MG approximation 
and to the resultant model which obeys the MG formula 
as stricto sensu (s.s.) MG model. We will denote all the 
quantities computed in the s.s. MG approximation with 
the script MG alone. 

When self-correlation is added to the s.s. MG model, 
deviations are obtained^ ' 10 ' 11 ' 32 . The resultant theory 
incorporates recurrent scattering. Differently to the MG 
model, even in the limit q£ — > no exact solution have 
been found and further approximations are needecU 2 -. 

We follow Ref. lj to derive the renormalized expressions 
of the single-particle polarizability, a, and the electri- 
cal susceptibility, x(q). The addition of radiative correc- 
tions amounts to the inclusion of self-polarization effects, 
closely related to the concept of local field. Because the 
self-polarization field depends itself on the dielectric con- 
stant, e e //, the problem becomes one of self-consistency. 
As a result, a is an implicit function of e e //- 

Because radiative corrections are affected by the 
spatial dispersion of x{q), further approximations are 
needed. Throughout this work we apply a quasicrys- 
talline approximation 20 . It has the convenient property 
that it is exact for the s.s. MG model in the effective 
theory limit. In particular, we will apply the s.s. overlap 
approximation 5 and a variation of it which incorporates 



recurrent scattering. 

We will restrict ourselves to purely radiative correc- 
tions. Non-radiative effects such as the collisional shift 
and Doppler and collisional broadening will not be con- 
sidered. Also, additional quantum correlations between 
the electrons of neighboring dipoles and the effect of 
the short-range interaction potentials between dipoles on 
their individual polarizabilities will be disregarde d 11 ' 25 . 

The paper is organized as follows. In Section [TT] we 
review our approach in Ref. [l| to derive an expression 
for the renormalized single-particle polarizability. The 
approximations used in the subsequent sections are de- 
scribed. In Section [IIII we apply the precedent formalism 
together with the overlap approximation to the compu- 
tation of the radiative corrections of the dielectric con- 
stant. In Section [IVI we compute the spectrum of coher- 
ent emission in the three dimensional s.s. MG model. In 
Section [V] we introduce recurrent scattering. We obtain 
corrections to the s.s. MG dielectric constant and to the 
spectrum of coherent emission. In Section |VT1 we propose 
both numerical and laboratory experiments where our re- 
sults can be tested. For the measurement of the coherent 
emission spectrum, formulae in two-dimensional samples 
are obtained. A numerical example is carried out. In the 
Discussion section we compare our approach with previ- 
ous works. Special attention is given to the role of the 
local field factors. 

About the notation, we label three-spatial-component 
vectors with arrows and three-by-three tensors with over- 
lines. We denote the Fourier-transform of functionals 
with g-dependent arguments instead of the r-dependent 
arguments of their position-space representation. 



II. BOTTOM-UP APPROACH 

Our procedure to compute the radiative corrections on 
the dielectric constant involves two complementary steps. 
In the first one, we renormalize the single-particle polar- 
izabilities. In the second one, we compute the electrical 
susceptibility. The former procedure gives rise to the so- 
called corrections over the Claussius-Mossoti formula 18 . 
The latter gives rise to the so-called corrections over the 
Lorentz-Lorenz or Maxwell-Garnett formula 3 ^. The fact 
that dipoles are treated as point-like allows us to per- 
form both renormalization processes separately. Self- 
consistency between both of them is required. 



A. Renormalization of the single particle porlarizability 

In a statistically homogeneous and isotropic medium, 
coherent emission of frequency lo is propagated from an 
external source by the Dyson propagator, G u (q). In 
terms of the susceptibility tensor, x"^), the transverse 
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and longitudinal components of G u '(q) read, 
Gl(q)= ' 

GTT(«) = 



k 2 [i + xl(q)} 
1 

fc 2 [l + X<f( g )] 



(1) 
(2) 



The propagator of the self-polarization field was found in 
Ref. Ill to be 



em 



^G ± (q) 
pa 



k 2 [l + Xl(q)]-Q 2 



pa 



Gii(q) = -r 



pa k 2 [l + X ^(q)] 



(3) 
(4) 



In the above formulae, it is implicit that x°i || (?) can 
be expanded as a series of one-particle-irreducible (1PI) 
multiple-scattering terms of order n, 



(«) = Ex£ii(« 



n=l 



E 

n=l 



W )p"fi», (5) 



where p is the average numerical volume density of scat- 
terers and the functions X^l(q, u) incorporate the spa- 
tial dispersion due to the spatial correlations within clus- 
ters of n scatterers. G u (r,r) is the propagator of virtual 
photons of frequency uj emitted and annihilated at the 
position f of a dipole. It amounts to the radiative correc- 
tions which give rise to the renormalized polarizability a. 
We define the 7-factors, 27" and 7^ as the traces of the 
transverse and longitudinal divergenceless components of 



2ll = 
1\\ = 



2d J fc 



Xl(q)/(pa) 



(2ir)nk*[l + xl(q)]-q- 
d 3 k r 1 Xf (g) 
(2tt)3 Lpafc 2 [l + xf(g)] 



(7) 



where G^\, (q, w) are the transverse and longitudinal 
components of the EM field propagator of frequency u in 
free space -i.e., Eqs. (jll2[) with x± 11 = 0. In the spatial 
representation, the longitudinal component corresponds 
to the electrostatic dipole field propagator, 



Gftat. 0, W) 



k 2 



(47rr) ! 



(8) 



while the transverse part corresponds to the radiation 
field propagator, 



G^ d .M 



A kr 



—Aur 



A kr 



1 



— Aur 



(9) 



Their respective traces contain divergences in the real 
axis which correspond to the terms subtracted in 
Eqs.jSIZl). In Ref. 1 we obtained in function of the 7- 
factors, 



a{k) = 



fc 2 aoK{7i 0) } + ^ 2 ao[2 7± +7||] 



(10) 



and parametrized a as a Lorentzian polarizability, 
a{k) = a' Q k 2 res [k 2 res -k 2 - irk 3 /(ck 2 res )}-\ 



(11) 



a are the renormalized values of the reso- 
nant frequency, the line-width and the electrostatic po- 
larizability respectively. Appropriate regularization of 
the divergences in terms of the bare resonant frequency, 
ujo = cko, and the bare electrostatic polarizability, ao, 
-3 

are, 



yields^ ^R.{2j^} = pr^. The renormalized parameters 



(o) 



2ir 

-r T2^{2 7 



(0) 



■2 7± 
-2 7J 



7||}|fc=*v, 
-7||}|fc=fc r 



and a' Q = a (k /k res ) 2 , 
where r = caokQ/6ir and k res satisfies, 

(k/k ) 2 - 1 = ^a k 2 $i{2 1± +7||}|fc=fc r 



(12) 



(13) 



(14) 



B. Correlation function and length-scales 

As mentioned in the Introduction, spatial dispersion ef- 
fects depend on the specific microscopical structure of the 
complex medium. While the dielectric constant e e ff is, 
by definition, the zero mode of the renormalized dielectric 
function e(q), the renormalization of the single-particle 
polarizabilities is affected by near-field effects where the 
microscopical details matter. As a first approximation, 
if the medium is not long-range-strongly-correlated, the 
Kirkwood superposition approximation 29 is applicable. 
That is, in good approximation all the 1PI diagrams can 
be computed out of the two-point correlation function. 
Generically, we can approximate 

h{r) ~ -f(r -0+ P^S^ir) + C^(r - £). (15) 

In this equation, f(r — £) accounts for the exclusion vol- 
ume around each dipole and its precise form depends on 
the interaction potential between pairs of scatterers. It 
tends to one for r < £ and to zero for r > £. Usual forms 
are those of a Lennard- Jones potential and a hard-sphere 
potential. With no much loss of generality here we will 
adopt the latter form, f(r — £) = 0(r — £). It behaves 
as a virtual cavity. The one-dimensional delta function 
in Ea. (fT5"|) stands for the overdensity of first neighbors 
around a given dipole. The constant C relates to the co- 
ordination number. It might be relevant for high-ordered 
media. The three-dimensional delta function stands for 
the self-correlation and it is ignored in the s.s. MG model. 

Let us set C = for the moment and study the rela- 
tionship between the length-scales of the problem. The 
purpose of this qualitative analysis is to make an appro- 
priate selection of the relevant diagrams in the cluster ex- 
pansion. The length-scales are^S, the wavelength of the 
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radiation, A; the dipole radius, a; the radius of the exclu- 
sion volume, £; the average distance between scatterers, 



C. The quasicrystalline and the overlap approximations 



P 



-1/3. 



and the effective length associated to the dipole 
strength, ~ jo^fe)! 1 / 3 . Off-resonance, the latter scale can 
be identified with that of the bare electrostatic polariz- 
ability, [ao(fc)] 1 / 3 . For 'classical dipoles' oto(k) is of the 
order of a 3 while for atomic dipoles, ao(fc) ~ k^ 2 r e , r e 
being the electron radius^. Because the dipole approx- 
imation demands £ > 3a and fc£ <C 1 for an MG di- 
electric, the neglect of a is generally justified under any 
circumstance. Hence, dipoles can be treated as Rayleigh 
scatterers. 

Off resonance, the scattering cross-section of single 
dipoles is much smaller than £ 2 . Recurrent scattering is 
weak and the most relevant correlation function is that 
of the exclusion volume, — 0(r — £). In contrast, at reso- 
nance the scattering cross-section of Rayleigh scatterers 
is of the order of A 2 es . Therefore, there is a change in 
the scale hierarchy when passing from the off-resonant to 
the resonant regime which affects the relation between 
£ and la(fc)! 1 / 3 . The influence of the exclusion volume 
becomes much weaker since the scatterers may overlap 
optically for |a(fc)| > £ 3 . The inherent self-correlation 
function p~ 1 S^ 1 \r) becomes relevant as it accounts for 
multiple-reflections between scatterers. That is, recur- 
rent scattering correction a 12 ' 32 ' 33 . 

In addition, another change in the scale hierarchy 
might appear at resonance depending on the relation be- 



tween X res and p~ ' . In the first place, if k res p~ /A < 1, 
the convergence of the series in Eq.© gets slower. Also, 
the dipoles overlap optically within a space where near 
field effects are dominant. Therefore, in good approxi- 
mation only the electrostatic dipole-dipole interaction of 
Eq.® is relevant for the recurrent scattering which en- 
ters the pair polarizability at resonance. This is con- 
firmed by numerical simulations 12,33 . In the opposite 
limit, k res p~ x ^ < 1, far-field interactions must be in- 
cluded based on numerics as well 3 ^. 

The above arguments and the numerical simulations 
of Ref. [ll|, [23|, and [32| make us to conclude that, leaving 
aside near-field corrections on the single particle polar- 
izabilities, the dielectric constant will be well approxi- 
mated by the MG formula off-resonance. Also based on 
the above arguments and the numerical simulations of 
Ref. [TH, 0, and Hfl strong deviations w.r.t. MG are 
expected at resonance. Note that the transition from 
the off-resonant regime to the resonant one takes place 
within a very narrow frequency range less than T. That 
is so since T -C uj res and fc res £ <C 1. Therefore, in good 
approximation, we can treat separately both regimes. 

In this paper we restrict ourselves to the off resonant 
regime and treat the recurrent scattering corrections as 
perturbations. We leave for a separate paper the compu- 
tation of the line-width and the resonant frequency shift 
of Eqs. (112ll4[) which require calculations at resonance. 
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FIG. 1. (a) Diagrammatic representation of the geometrical 
series which amount to Xmg at q = 0. (b) Series of 1PI dia- 
grams which amount to x' 3 ' in the Kirkwood approximation. 
Only the exclusion volume correlation function enters the MG 
model, (c) In the long- wavelength limit, — > 0, the series 
in (b) is equivalent to the diagram which contains consecutive 
correlations onl y 31 ' 34 . For shorter wavelengths, (6) and (c) 
are related via the overlap approximation^. 

It was proved in Ref. |8| and |3fJ that, in the long- 
wavelength limit, g£ — > 0, if only the negative correla- 
tion of Eq. (fT5l) is considered, h MG (r) = —f(r — £), the 
effective susceptibility is exactly that of the MG formula 
provided that /(r — £) is isotropic^, 



pa 



(X 



MG 



3) 



(16) 



The 1PI diagrams which amount to x(g) are of the kind 
of those for \^ (<?) depicted in FigJTJ 6) in the Kirkwood 
approximation. In the limits g£ — > 0, fc£ — ¥ and for 
h M (r) = —f(r — £), the series is equivalent to that in 
FigJIJa) which contains only consecutive correlations. 

Generically, if that equivalence is extended for any 
h(r) to shorter wavelengths, the relation becomes an ap- 
proximation which is referred to as quasicrystalline ap- 
proximation 20 (qc). In this approximation, the series of 
x(<7, w) becomes geometrical and the only quantity to be 
computed is x^ 2 )(g,o;), 



pa 



1 - Xx\(Q,v)/(p6i) 



(17) 



For the computation of the dielectric constant, only the 



zero-mode matters, xTffi 10 ) 



However, spa- 



tial dispersion is unavoidable in considering radiative cor- 
rections since shorter wavelengths of the above formula 
are needed for the computation of Eqs. (|6l7l) . 
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In Ref.l8landl34l the authors worked over a hard-sphere 
model in the quasicrystalline-MG approximation with 
h ov (r) = — 0(r — £). They referred to the resultant ap- 
proximation as overlap approximation -see Figfljc). In 
Section IIIII we apply this approximation to the compu- 
tation of radiative corrections. Because self-consistency 
is demanded there, we will refer to it as self-consistent 
overlap MG approximation. We will denote the quanti- 
ties computed in this approximation with a script MG " 
and we will use the script MG alone when the computa- 
tion is independent of the precise form of f(r — £). Since 
it is for q£ — > that the overlap MG approximation be- 
comes exact, the formula 



Xmg 



pa 



i-x 



(2) 

±,\\MG 



, v (q = 0,u)/(pa) 



(18) 



is exact and equals Eg. (1161) . 

Later on, we add the self-correlation term to h ov (r), 
Kec.( r ) = -®( r - + P^S^Hr). Its physical effect 
is the inclusion of corrections due to recurrent scatter- 
ing (rec). In this case, the quasicrystalline approxima- 
tion is not even exact in the long- wavelength limit. The 
reason being that the recurrent scattering diagrams in- 
volve multiply entangled integrals. Nevertheless, we stick 
to the quasicrystalline approximation in Section [V] Be- 
cause still we take f(r — £) = 9(r — £), we will refer 
to the resultant approximation as self-consistent overlap 
recurrent approximation. We will denote the quantities 
computed in this approximation with a tilde at the top 
-except for a in which the tilde denotes the renormalized 
value of a. 



III. COMPUTATION OF e e// IN THE 
SELF-CONSISTENT OVERLAP MG APPROXIMATION 

As mentioned above, in the s.s. MG model the only 
relevant correlation function is that of the exclusion vol- 



ume. In such a case Eg. ([16 
with Eq.([TU]), we obtain, 



is exact and, in combination 



pa 



'MG 



1 - 



\pa 



= 1 + pa k 2 



1 



' Xmg 



1 



(19) 



+ -^ a o k o k ( 2 7i + 7||)fc ~ ^ aokoP 

In the following, we calculate consistently the 7-factors in 
the overlap approximation. The resultant formulae are 
function of £ and xmg- Those formulae together with 
Eg. (1191) give rise to an implicit function of the complex 
variable xmg- 

Let us compute first 2~/± . Applying Eq. @ with x± (l) 
given by Eg. (|17p . the transverse polarization propagator 
reads, 



G?>(«) 



1 



-pa[k 2 G { l\q) 



(20) 



The corresponding 7-factor is 
r C 2 



2 7 r 



c 2 - 



(21) 



n=0 



(^)- 2 x ( L 2) (C,Q)] n }-5R{2 7 [ 0) }, 



L C 2 -Q 2 

where Q = qt;. Hereafter we omit the explicit depen- 
dence on uj — ck unless necessary. Using h ov (r) for the 
(2) 

calculation of x\ (C> Q) we get, 
(ap) 2 k 2 



X±m G ov (q) = 



d 6 re lq r Q{r - £) 



x Ty{G (0) (r)[I - q ® q}} 



(ap) 2 k 2 f dV 



2 J (2tt) 3 
+ G ( L °V)cos 2 + G\°\q') sin 



h° v (\<r-<j\)\Gf\q') 



(22) 



where q is a unitary vector parallel to q and h ov (Q) — 
-4< 3 ji(Q)/Q = 4< 3 Qcos ^ sinQ with j 1 the spherical 
Bessel function of first order. At leading order in £ = fc£, 



X±mg°- (Q) - (P & ) 



> h(Q) 
Q 



o(( 2 ). 



(23) 



Inserting Eq. ([23|) in Eq. (|2T|) . we can decompose the re- 
sultant formula into long and short wavelength terms up 
toO(C°), 



27. 



MG° V 



Z 1±MG "T Z 7_LJ\/G°" 1 



7 00 -oo 

^i F MG = 4- c {^-P^ 1 d Q] 
' ^ n=0 J ° 

L/3]"}, 



with 

Q 2 



L ( 2 -Q 2 



2 7±MG 



u 00 r°° 

n c Jo 



(24) 



(25) 



Eq. ([24|) carries the far field (FF) propagating modes 
while Eq.(|25|) contains the near field (NF) contribution. 
By taking the limit £ <C 1, it is immediate to identify 

in Eg. (|24| the term —1/3 (pa) 2 with the zero mode of 

(2) 

~X±MG ov (Qy W e take advantage of this to arrange the 
0(C°) leading terms of Eq. ([24"l) in an exact expression, 

1 ff a \ k e M G + 2 . , . 

2 l±MG[ k ) = q Vcmg- (26) 

Z7T O 

27jj^ G o« amounts to those non-propagating short 
wavelength modes which scale as l/£. The corresponding 
series is rapidly convergent. We give below the first five 
terms of that series, 



Z I±MG 



-k 
2^C 
334 
31185 



l -pa-^{pa) 



47 

+ 7^7>«) 3 (27) 



1280 



(pa) 4 



6891623 
2145927168 



(paf + 
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Let us consider next the computation of 1^ IG ■ The 
propagator of longitudinal modes in the quasicrystalline 
approximation reads, 



give a closed expression for Q 



MG" 



(Q), 



-,MG m ' 
'II 



{ £mg - 1 hXQ) 

'ea/g + 2 Q 



/(C)]- 1 . (31) 



Gl c (q) = ^[l + P&-xf\<l)/p&} 



In the self-consistent overlap approximation, 



(28) The use of this formula in the calculation of 7j| im- 
plies the computation of an infinite series of roots of 
1 - 3 eMG 7o J1 ^ Q) /(C) 8 - Instead of that, we rather ex- 
pand Eq. (PT|) in powers of pa as this allows us to keep 
control over the order of accuracy of the series, 



xfL°M = (Pa) 2 k 2 j d 3 r e ^ ? 6(r-0Tr{GW(r)[g®g]} r MG°" 

d 3 q' 



-(pafk 1 



:h ov (W-q\) 



(2tt) 3 

g(,°V)cos 2 + G < f ) (q') sin 2 



The above equation can be computed in closed form, 

MQ) 



Finally, we arrive at 
(29) > i oo 

E(- 1 ) n+1 (^)"/ n (o 



(32) 



MG ov 



2tt 2 C 3 

dQ Q 2 MQ)IQ] n . 



n=l 
l2r 



(33) 



Q 



-/(C)], /(0=2*e«(<+C). 



XiiM G ^(g) = (p«) 2 [i+ 

(30) 

Inserting Eq. ([50|) in Eq.J^EJ) and using the MG relation- 
ship for £mg m function of pa -Eq. (|18[ ). it is possible to 

I 



Note the similarity between the above equation and that 
for 2 r y I [ MGov in Eq. (|25|) . 7^ fG contains however addi- 
tional terms of order 1/C 3 and order zero in £ as a result 
°f /(C) m Eq. (|501) . We give below the first five terms of 
the orders 0(C~ 3 ), 0(C~ 4 ) and O(( ), 



7| 



MG" 



l ^[\pa + \{paY + ^{pa)U 



272 
2835 

544 
2835 



(pa) 4 
(pa) 4 



1088 
8505 



(pa) 4 



40949 
870912 
204745 



1741824 
204745 



(pa) 5 4 
(pa) 5 



2612736 



(pay 



(34) 
(35) 
(36) 



Because 2j^ MG contains only long wavelength modes, 
the overlap approximation is there as good as for the 
computation of xmg a t order 0(C°)- Hence, the super- 
script MG instead of MG 011 in Ed.([2l)]). On the contrary, 



I 

-f MG ° v together with 2^ MGov contain the contribution 
of short wavelength photons which probe the microscopic 
structure of the dielectric. The application of the overlap 
approximation in them is only exact for the O(pa) term. 
In function of Xmg-, the 7-factors read, 



2j MG- 



MG" 



. k e M G + 2 



2tt 
ik .1 

2^h 

--I 

k 1 

[Xmg 



Xmg 



Xmg 



V WG 

■X 2 mg ~ 0.051 X 3 MG + 0.055 X 1 /G - 0.015 X mg + ■••] 



Xmg + 0.105 X 3 MG - 0.027 X t /G + 0.006 X mg 



;Xmg 



; xif G +0.07xt/ G -0.03 X 5 A fG ...]. 



(37) 
(38) 
(39) 
(40) 



I 

As advanced, the 7-factors in Eq. ([57)l are function of £ in Eq. (fT9"|) . the computation of Xmg becomes a problem 
and xmg- Because xmg is also a function of the 7-factors 
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of self-consistency. 



V. BEYOND THE MG MODEL 



IV. THE SPECTRUM OF MG POLARITONS AND 
COHERENT EMISSION 

It was suggested in Ref. 0] that the coherent radiation 
emitted by one of the dipoles in a complex medium could 
be collected in the far field. According to the standard 
Beer-Lambert law, 



W 



Coh. 



r '-f\) = M/ Co "-(f)exp[-2fc|r' - r\n], (41) 



where k is the extinction coefficient, n — S{^/6mg} and 
W (r) is the coherent power radiated from the emitter 
site, f. For a statistically homogenous medium, 



W 



Coh. 



(f) = W 



d <! 25R{^4^}3{Gi(<z)}, (42) 



(27T) 



pa 



where W a = g^|pb| 2 , being p = e Q aE$(r) the dipole 
moment induced on the emitter by an external field of 
frequency w, Eq (r). The integral of Eq. (|4"2l is, up to 
prefactors, the spectral function of transverse polaritons, 

7rc J (ziry pa 

(43) 

The poles of 3{Gj_(g)}, kj° 1 ', satisfy the dispersion rela- 
tions for transverse normal modes, k 2 e±(q) — q 2 \ q=k Poi. — 

0. The residue ^{ x±{ ^'"' ) } was interpreted in Ref. 1 
as a renormalization factor. For an MG dielectric, the 
condition fc£ <C 1 precludes the existence of geomet- 
rical resonances and G±(q) contains a simple pole at 

k{\ oL = y/ma. In the s.s. MG model, jft{ 2^M!!l } = 
?R{eMG}+2 ^ yjYnch is a Lorentz-Lorenz local field factor. 
Thus, Eq.(03l) reads, 



LDOS 



Pols. 
MGX 



M = 



oj 2 K{e MG } + 2 



4ir z c 



2„2 



^{V^}- (44) 



We do not address here the longitudinal polaritons as 
they require further discussion^. They are irrelevant in 
the present context. 

Making use of the MG formula we can write po as 
an effective spherical dipole of volume p~ l and dielec- 
tric constant 6mg, Po = ^oP^ 1 f^j| -^o (f)- Thus, we 
obtain, 



K{e MG } + 2 



(45) 

l27 Jr?r \E%{r)\ 2 . Note that Eqs. (|44l45j) arc 
not subject to the overlap approximation but just to the 
MG approximation. 



where W,°, = 



In this Section we correct the MG and MG"" formulae 
with recurrent scattering terms due to the addition of 
the self-correlation which is ignored in the overlap MG 
approximation. 



A. Off-resonant recurrent scattering 

The overlap recurrent approximation consists of using 

^rec( r ) = ~~ ®( r ~ + P -1 <^ 3 H^) an( i applying the qua- 
sicrystalline approximation. The variation of the effective 
MG susceptibility due to the introduction of recurrent 
scattering in the long-wavelength limit has been stud- 
ied by a number of authors^ ' 10 ' 11 ^ 32 . We include short- 
wavelength variations. It is worth noting here a subtle 
point. Even though the recurrent scattering processes 
(denoted with the script rec.) are formally equivalent to 
some of the virtual self-polarization diagrams which cor- 
rect the bare polarizability, the physical content is differ- 
ent. That is so because the recurrent processes amount 
to the multiple reflections that actual probe photons, in- 
stead of virtual, experience within clusters of scatterers. 

The two-scatterer diagrams of XsUsc. are depicted in 
Figj^l They include an infinite number of reflections be- 
tween two particles. As it was explained in Section IlIBl 
we will consider separately the purely electrostatic (near 
field) contribution and the purely far field (scalar-like) 
contribution. In these two limits, we have the respective 
simplified propagators, 



G^(r) = (l-3ff)/47rfcV, 
G { °2(r) = [-e ifer /47rr]I. 



(46) 
(47) 



The corresponding 'two-scatter recurrent susceptibilities' 
read, 



= (2)™.,-, s_ (p^f MO) 
Xst.,sc. V^) ~ TJ—^st.^c.V) 



(48) 



^(k 2 a) 2m [G { 2, sc .(^] 2m [l-e(r-0l 



For the electrostatic contribution, 
{pa) 2 



Xlt. \r,u)) 



4?r 



■[l-9(r-0] : 



-(a/An) 2 



3ff 



r 3 (r 6 - (a/4n) 2 ) 
3r 6 - 4(5/4tt) 2 



r 6 - 4(a/47r) 5 



(49) 



The zero-mode of the above formula yields a correction 
over the factor l/3(pa) 2 of the usual MG formula. At 
leading order, 



x { ?^::(Q = o)^-(p&) 2 (&/^e) 



-3\2 



(50) 
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As expected, the longitudinal and transverse zero modes 
are equal. Beyond the long-wavelength limit, 



x^r(Q)-(p«) 2 #(-/4^ 3 ) 2 , 



(51) 



~<2)rec. 

x\Jt. 



(Q) ~ -2 X fl r t ec -(Q) + ( P ay — (a/4ner (52) 
cos Q(120 - 6Q 2 + Q 4 ) + (720 
24Q 2 + 2Q 4 - Q 6 ) + Q 6 CosIntegralQ 



We find that, off-resonance, the dimensionless expansion 
parameter for the electrostatic recurrent susceptibility 
is a /An!; 3 . This was already noticed by Kirkwood and 
YvoniiiS. The precise computation depends on the type 
of interaction between the scatterers. For a hard-sphere 
potential Alder et al^ and Cichocki and Felderhof 32 
computed numerically xlt c ' (Q = 0) as a series expansion 
on z = 2 3 a/47r^ 3 . They found for the 0(z 2 ) term a nu- 
merical factor in the range [0.016,0.028] which depends 
slightly on the filling factor. Our analytical calculation 
yields 1/2 6 ~ 0.016. 

For the far-field scalar-like contribution, 



Xi 2 c )rec -(r» 



-(pa) 2 k 6 a 2 [l-e(r-0] 



l-fc*d 2 [GS2(r)] a 
Its zero-mode at leading order is, 



(53) 



xi 2 c )rec -(Q-o) 



- (pa) 2 k 6 a 2 /4ir [ lEu - 1 + lg 3( - »tt/2] . 

(54) 



Beyond the long- wavelength limit, 

xi 2 c )rec '(Q) ^ Up & ) 2 k 6 a 2 /4n\2(l - lEu ) ~ \g(% 2 Q 2 ) 



2 

3C, , K-Q , , . , QA/n 
Q lg( 3CTQ ) + i7r(3C/Q 



|i - 3C/QI) 



where 7£ U is the Euler constant. We find that the dimen- 
sionless expansion parameter for the far-field recurrent 
susceptibility is ak 3 . This parameter is much smaller 
than a/47r£ 3 in an MG dielectric. Therefore, the discard 
of the scalar contribution is justified in good approxima- 
tion. 

I 



■pa 



(2)r« t _ I I I/—-., 1 I + 

1 i-M M M it 



FIG. 2. Diagrammatic representation of the two-scatterer re- 
current susceptibilities of Eq.l 



B. Computation of e e // in the self-consistent overlap 
recurrent approximation 

We proceed to compute the modified 7-factors. To 
do so, we include the recurrent terms of Eqs. (|50H52"j) in 

x {2) (Q^), 



* (2) (Q, <*>) = xfloo, (Q, «) + x^ rec - (Q, w) 



-,(2 



~(2)rec. 



(55) 



and apply the quasicrystalline approximation which leads 
to x(Q,w) via Eq. ((T7| . By comparing Eqs. (150H52j) with 
those corresponding to the overlap MG approximation, 
the computation of X (Q^) just involves the substi- 
tution of the factor (pa) 2 by (p<5) 2 [l + (a/47r£ 3 ) 2 ] in 
(2) 

Xmg ov (Q' IjJ )- ^he resultant formulae are exact, in our 
approximation, up to 0(a 4 ) in Xe ff and up to 0(a 3 ) in 
the 7-factors. Higher order accuracy can be obtained if 
needed out of Eq. (|4"5|) . The relation between pa and X ef f 
gets modified as 



Xeff = pa 
or equivalently, 

3 Xe// 



3\2i 



1 - -pa[l + (a/M 6 ) 



(56) 



pa : 



Xeff + 3 



1 - 3(4^ 3 P)" 2 X 3 



eff 



x [9(4< 3 p)- 2 x 3 f f + (xeff + Sf/Sr 1 . (57) 



Using the above relations, the modified 7-factors read, 



2tT(*) 

27 r 



k e 



eff 



2r 



2tt 3 

-k rl - 



2ttC 

-k 1 r 
k 1 rl 

~ 2^C 
k 

— i 



1 - 3(4^ 3 p)- 2 Xe//[9(4^ 3 p)" 2 Xe// + (Xeff + S) 3 /^^ V^ff , 

1 - . >_0 Q ,14 

15 



3 _ 2 



(47re 3 p)- 2 + 0.181]x 3 



■(4 7 re 3 /9 )- 2 + 0.11] X 4 // . 



2tt 



xeff + \xlff + K^epy 2 - \\xlff + omzxiff 

Xeff + \xlff + [^(4< 3 p)- 2 - 0.076JX 3 // 

1 



Xeff + 3X 2 // + l^< 3 Pr 2 - 0.05JX 3 // 



(58) 
(59) 
(60) 
(61) 
(62) 
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There is just a subtle point about the 0(£ x ) and C(C°) 
terms of order (a/4ir£, 3 p) 2 x 3 eff in Eqs. ([59l61l62p . Differ- 

(2) _____ ______ 

ently to the expression for x\mqov m Eq. (j3"0")) . Eq. (j5"2")) 
does not contain terms of orders 0(£ 2 ,£ 3 ). Therefore, 

(2) (2.) 

the substitution of Xn (Q) instead of X\\mg°v(Q) m 
Eg. ([28)) does not yield directly all the terms of orders 
O(C _1 ,C ) in 7||. The reason being that in the formu- 
lae for Xsi.sc? we have disregarded the coupling between 
electrostatic and radiative bare propagators. However, 
exploiting the similarity pointed out after Eq. (|33[) be- 
tween, in this case, 711 and 2 ; y^ F , those terms can be 
derived. The coincidence between the longitudinal and 
transverse terms at order pa(a / 47r£ 3 ) 2 (equivalently, at 
order (a/47r£ 3 p) 2 Xe//) is depicted diagrammatically in 
FigEl 

It is worth noting that, for a given order in a or 



\=-0K) 



T(r) 



'■ pa 



FIG. 3. Diagrammatic representation of the equivalence be- 
tween the transverse 7-terms (l.h.s. diagram) and the longi- 
tudinal 7-term (r.h.s.) of orders 0(C _1 ) and C(C°). 

Xeff, those terms which include recurrent scattering fac- 
tors a/47r£ 3 may be greater than those involving just pa 
factors. The weight of a/47r£ 3 relative to pa can be com- 
puted as a function of the effective filling factor, 0. Treat- 
ing the scatterers as hard spheres of radius £/2, regardless 
of the internal radius a of the scatterers, p — Thus 
we find 4ir£ 3 p = 240, which is less than one for <fi < 0.05. 



C. Correction to the MG coherent emission 



Recurrent scattering introduces also a modification 
on the local field factor w.r.t. that of Lorentz-Lorenz. 
Correspondingly, the expressions for LDOS^ ' s (cj) and 
W Coh -(f) get also modified w.r.t. Eqs. (|44l45j) . Using the 
relation of Eq. (|57|) . they read now, 



LDOS 



Pols. . . 
MG±\ UJ ) 



4ir 2 c 2 



5R- 



i^lL±l\l 3(4 7 reV)- 2 ^ // [9(4< 3 p)- 2 ^ // + (Xeff + V 3 /^ 1 } "^{v^}, (63) 



W 



Coh. 



eff 



e eff 



1 - 3(4 7 re 3 p)- 2 ^ // [9(4< 3 p)" 2 Xe// + (Xeff + 3) 3 /3]- 1 



x 5R{f££L±^ [1 - 3(4^ 3 p)- 2 x 3 // [9(4^ 3 p)- 2 ^ // + (Xeff + V 3 /^ 1 ] ^{V^Ffh 



(64) 



VI. EXPERIMENTAL PROPOSALS 

In principle, numerical simulations based on the 
Coupled-Dipole-Method (CDM)! 6 ^- or the Boundary- 
Element-Method (BEM)^1 could test our analytical re- 
sults. However, to the author's knowledge, there is no 
numerical code that takes into account consistently all 
the radiative corrections on a and Xeff- Some numer- 
ical works consider only electrostatic corrections due to 
actual recurrent scatterin g 12 ! 26 . Others do include ra- 
diative corrections but only those in free-space 23 . Some 
authors focus on deviations w.r.t. the Claussius-Mossotti 
formula due to the finite size of the dipole constituents 3 —. 
An improvement was made by the authors of Ref. [TH 
who carried out an a priori renormalization of the po- 
larizability of point dipoles in a lattice by including the 



lattice dispersion relation. Further on, the authors of 
Ref. [3^ have proposed a modification over the original 
CDM formulation which accounts for the integration in 
a of the near-field corrections within a lattice cell. Other 
works include additional radiative corrections on x&ff- 
In Ref. |34| retardation effects for ( ~ 1 are considered. 
In Ref. [2J the corrections come from the finiteness of the 
macroscopic sample. 

In the following, we propose both numerical and labo- 
ratory experiments to test our results. In the first place, 
we intend to verify the accuracy of the self-consistent 
overlap MG and recurrent approximations in the com- 
putation of the radiative corrections on a sample made 
of hard-sphere scatterers. Second, we aim to verify the 
generic form of our formulae for Wg° G l -{r) and W Coh {r) 
in an MG dielectric. The latter formulae are the most 
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universal result of this article as they are nearly indepen- 
dent of the microscopical details of the medium. That 
is, they are the result of considering only the inherent 
correlations -i.e., exclusion volume and self-correlation- 
whose detailed form is almost irrelevant for long wave- 
length propagating modes. 



A. Test of the approximations 

1. Numerical experiments on emission 

Despite the fact that numerical simulations do not 
take proper account of all the radiative corrections, the 
formula derived in Ref. [H for the total power emit- 
ted/absorbed by an induced dipole is still testable by 
means of CDM numerical experiments. It was found in 
Ref. [H that 



W, 



Tot 



^|a^| 2 [3{2 7 f+2 7±+7|| } (65) 



3 3{a } 



k 2 



If radiative corrections were properly incorporated in the 
renormalized polarizability, the above equation would be 
equal to ^\Eq | 2 3{<S}. Such an equivalence is a con- 
sequence of the optical theorem in a complex medium. 
Nonetheless, even if radiative corrections are not consis- 
tently incorporated in a, Eq. (|65[) can be tested numer- 
ically in aggregates of point dipoles like those of Mal- 
let, Guerin and Sentenac 23 . The reason is as follows. 
In Eq. ([3T|) and Eq. ([S"5f the 7-factors can be written in 
function of xmg and Xeff because Xmg and Xeff are re- 
lated to a by Eq. ([TBI) and Eq. (|57|) respectively. Likewise, 
we can use Eqs. (|16l57[) to replace cxEq in Eq. (pSl with a 
function of xmg > Eq (r) or Xef f > Eq (r) respectively. That 
is irrespective of the manner a is renormalized. In par- 
ticular, a can be taken to be the in-free-space value, a, 
as used in some simulations 23,28 . By exciting one of the 
dipoles in the aggregate with an external probe field as in 
Ref. and HH, integration of the total intensity in the far 
field must adjust to the first term of Eq. (p5|) . with a re- 
placed with a and disregarding microscopical absorbtion 
in the medium. The latter means that the macroscopic 
absorbtion is just due to microscopic incoherent scatter- 
ing in agreement with the interpretation of the authors of 
Ref. HD- As mentioned above, the ^-dependent terms are 
only exact at 0(1) in XMG, Xeff ■ Therefore, this would 
test the accuracy of the self-consistent overlap MG and 
recurrent approximations in the near field. 



2. Laboratory experiments on transmission 

Direct verification of our approach needs of laboratory 
experiments. Accurate information of the microscopic 
parameters a®, ko and £ is required as well as isolation of 



the radiative effects from those which are non-radiative. 

Experimental measurements are performed on the 
propagation of a normally incident light beam through 
a dielectric slab. The refractive index, n, and the ex- 
tinction coefficient, k, are derived out of the formulae of 
the transmition/reflection coefficients. Because the ef- 
fects of radiative corrections are expected to be more 
relevant under quasi- resonant conditions, instead of ab- 
solute measurements of n and k experiments concentrate 
on their spectra close to the resonance. The signatures 
of the radiative effects under investigation are the reso- 
nant frequency shift and the band broadening w.r.t. the 
values of the isolated dipoles in free space. According 
to our formulation, the frequency shift must be given 
by Eq. (fT3]) plus the Lorentz-shift, — ^paok 2 , in the MG 
model. The latter must be corrected by an additional 
term — ^paok 2 (j^j) 2 when recurrent scattering is con- 
sidered. The line- width must adjust to Eq.(fT2"j). 

A pioneering experiment on this issue is that of Ref. [4^ 
on an atomic gas of potassium. There, selective-reflection 
techniques were employed to measure the frequency shift 
of the transition 4 2 Si/ 2 5=^ 4 2 Pi/2- The problem to 
test our formalism in that setup is that both the Doppler 
broadening and other collisional effects dominate over the 
radiative effects. A good estimate for £ is also difficult to 
obtain. For cold atoms it would be of the order of twice 
the van-der-Waals radius. At high temperature however 
it is affected by the collision cross-section of the atoms. 

On the opposite extreme is an atomic glass. It has 
the advantage that £ can be accurately determined and 
vibrational effects well controlled. However, atomic or- 
bitals overlap in such a way that the electronic band con- 
figurations can be very different to those of the individual 
atoms isolated. That makes the estimate of ao difficult 
to achieve as its value may be very different to that of 
isolated dipoles. On top of that, there might be contri- 
butions of free electrons, which introduce an additional 
source of uncertainty in computing the component of the 
susceptibility due to plasma oscillations^. 

We conclude that suitable experimental setups must be 
such that undesired collisional effects reduce to a mini- 
mum; control must be kept over the parameter £; and 
the electronic structure of the dipole-constituents must 
remain unaltered either by their mutual cohesive inter- 
actions or by their interaction with a substrate. All these 
conditions meet in 'artificial solids' like those prepared by 
Liz-Marzan, Giersig, Mulvaney and Ung^r— . There, the 
dipole constituents are gold nanoparticles. The particles 
are first coated by silica shells prior to their layer-by- 
layer assembly. That way, the nanoparticles do not form 
aggregates which could modify their single-particle elec- 
tronic structure otherwise. Because the metal is nobel, 
no chemical reaction exists with the substrate. As a re- 
sult, the dielectric properties of the isolated gold particles 
are well described by Drude's model prior to formation 
of the solid and a good estimation for £ can be obtained. 
In Ref. |44] the authors carried out experimental measure- 
ments of the absorbance and the plasmon resonant shift 
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in thin films as a function of the filling factor. The ex- 
clusion volume is well approximated by twice the silica 
shell-thickness plus the radius of the Au particles. The 
density and size of the particles are known by construc- 
tion. The results are compared with the ordinary MG 
formulae. The agreement is qualitative but not quanti- 
tative enough. Since our approach incorporates the near 
field terms which are missing in the ordinary MG for- 
mula and the coated particles are perfectly described as 
hard-spheres, we expect to obtain a better fit. However, 
because the present work covers only the off-resonance 
part of the spectrum of n(ui) and the at-resonance 

curves of Ref. [44] cannot be compared with our formulae. 
Nevertheless, a similar experimental setup might still be 
useful to verify our formulae for the spectrum of coherent 
emission as explained below. 

B. Test of the Coherent Emission Spectrum 

For practical purposes, the measurement of the coher- 
ent intensity is of interest in signal processing and trans- 
mission at both macroscopical and microscopical scales^. 
Also, there is a fundamental reason why the verifica- 
tion of our result is of interest. That is, the formulae 
of W^g '(r) and W Coh -(f) are nearly independent of the 
microscopical details of the medium for low frequency 
emission. Therefore, they can be considered as generic 
for the kind of media which adjust to the so-called virtual 
cavity scenario. Those are, media made of a collection of 
equivalent dipoles in which the emitter is itself one of the 
host dipoles. The verification of our result would throw 
light into the historical problem of the role of the local 
field factors in dipole emission. 

1. Numerical experiments 

The formula of Eq.([i5|) for W^a^r) is exact for the 
s.s. MG model since the overlap approximation is exact 
in the long-wavelength limit, — > 0&2£i. When correc- 
tions due to the inherent recurrent scattering are incor- 
porated, the formula in Eq. ([M|) for W Coh (r) is exact 
up to order 0(xl/f) f° r a hard-sphere model. Because 
both W^jQ'ir) and W Coh -(f) depend only on macro- 
scopic quantities, Eqs. (|45l64[) do not need of the knowl- 
edge of any microscopical parameter to be validated^. 
In numerical simulations it is obvious that the access to 
microscopical parameters is not a problem. On the con- 
trary, in experimental setups it is, and that is why the 
validation of Eqs. (|45l64[) would be easier than the direct 
validation of tuc and €e// ■ We sketch first how to ver- 
ify Eqs. (|45l64p in the numerical simulations with dipole 



a Except for the dependence of W Coh '(r) on £ which can be re- 
placed by its relation with <j>. 



aggregates of the kind of those in Ref. [H and [H. As for 

Tot _ 

, the dipole at r must be excited by an external field. 
Then, it is necessary to compute statistically the average 
coherent field. Performing the statistical average at a 
single point should be enough provided that the orienta- 
tion of the emitter is randomized as well^S. Differently to 
the measurement of the total radiation which can be per- 
formed out of a finite size sample, the point f 1 where the 
coherent field is measured must lie well inside the sample 
and in the far field with \r' — t\ 3> £■ Randomization 
of the emitter orientation implies that the intensity as- 
sociated to the average field at f* is the angular-average 
intensity at a distance \r' — r\ from the emitter. Thus, the 
square of that field times 27r|r' — r\ 2 ceon is the quantity 
to compare with the Beer-Lambert formula of Eq. (|41[) . 

2. Laboratory experiments 

The preparation of a suitable laboratory experiment 
to test Eqs. (|45l64[) involves a number of difficulties. 
In principle, independent measurements of the refrac- 
tive index, the extinction coefficient and W Coh (r) could 
be performed. To determine n, K, standard transmit- 
tion/reflection measurements on monochromatic light 
can be carried out. For the computation of W Coh ^r), 
a unique dipole must be excited and the far field coher- 
ent radiation of Eq. (|4T|) must be detected. For the latter, 
the detector must be introduced in a manner that it does 
not modify the surrounding dielectric medium. For in- 
stance, it can be placed right at the edge of the sample. 
Coherent radiation can be filtered using interfcromctry 
techniques. Unfortunately, there exists a potential diffi- 
culty in the excitation of the emitter. Should the emitter 
be an isolated particle in free space or a single particle 
doped with a fluorescent ion, a laser beam could be used 
as in Ref.H3. However, the stimulation of a unique dipole 
in an MG dielectric using an incident light beam is not 
possible due to the inherent diffraction limit. That is, a 
monochromatic light beam cannot resolve distances less 
than ~ A/ 2. Therefore, it would be necessary that the 
radius of the exclusion volume, £, be at least of the or- 
der of A to avoid the excitation of more than one dipole. 
However, the low frequency requirement implies ^/Acl 
and so both conditions cannot be satisfied at the same 
time. As an alternative, we propose the use of near-field 
optical techniques^. That implies the introduction of an 
external device inside the dielectric, namely the metallic 
tip of a near- field-scanning-optical- microscope (NOMS). 
The tip excites the emitter dipole via near-field induc- 
tion. In a three-dimensional dielectric that would break 
manifestly translation invariance since the tip would act 
as an impurity. That is, the tip could scatter part of 
the emitted radiation and also interact with the dipoles 
around the emitter. Thus, unpredictable deviations from 
Eqs. (|45l64p are expected. To go around this problem, we 
rather propose to work on a quasi-two-dimensional di- 
electric. The tip must be placed perpendicularly above 
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the sample as in FigfJ] so that it does not perturb the 
medium. In order to avoid the interaction of the emitted 
radiation with the tip, only on-plane polarized modes - 
i.e., p-polarization or TE-modes- are of our interest. An 
analogous formula to that in Eq. (|45|) is worked out in Ap- 
pendix |X] for the TE coherent emission per unit length, 
I, 



NSOM 



-MG 



KM/ G } + 1 



(66) 



7r 1 arctan 




Far-field detector of coherent emission 

FIG. 4. Schematic representation of the experimental setup 
suggested in Section TVI B 21 The p-polarized near-field of an 
NOMS is used to excite one of the gold nanorods of radius a 
and length I, with k' 1 > ! > a . The on-plane polarized ra- 
diation (TE-modes) experiences multiple-scattering processes 
until collected in the far-field. The detector is placed right at 
he edge of the sample at a distance r' away from the emitter. 
The coherent signal is filtered using interferometry. 

l^Wp h '{f) can also be found which captures the ef- 
fects of recurrent scattering up to 0[{x P e j f) 5 ]- We pro- 
ceed in the same manner as for W Coh \r) in Eq. (f6"4")) . The 
two-dimensional expressions analogous to Eqs. (|56l57l) are 
respectively, 



Xeff 



pa 



2\2l 



1 - ^cra p [l + (dp/27r£ ) 



(67) 



where the script p stands for p-polarization and = 
|£^(r)| 2 with cr being the numerical surface den- 



16c 2 CT 5 

sity of dipoles. The preparation of two dimensional sam- 
ples could be achieved using the same coating techniques 
as for three dimensional artificial solids^. Alternatively, 
regular arrays of the kind of those in Ref. can be used. 
Beyond the MG approximation a formula for 



2X P 



eff 



x [6(2< 2 CT )- 2 (x: // ) 3 + {xl ff + 2) 3 /2]- 1 ] • (68) 

And so, by replacing the local field factor with its 
recurrent-corrected value we end up with, 



%ff + 1 



1 - 2(2<V)- 2 (^ // ) 3 [6(2<V)- 2 (^ // ) 3 + {Xlff + 2) 3 /2]- 1 



x 5R{^^ [l - 2(2^ 2 a)- 2 (^ // ) 3 [6(2^)- 2 (x P ef f) 3 + (Zff + 2) 3 /2]- 1 ] 



1 + 7r arctan I 



(69) 



The verification of Eqs. ()66|69p in numerical simulations is 
also computationally less demanding than those in three 
dimensions -see eg. Ref. Hi] and [H. 



C. A numerical example 



We perform here the numerical computation of the 
dielectric constant of an MG dielectric made of classi- 
cal hard-sphere dipoles. Also, we compute the power 
emitted by one of the dipoles when stimulated by an ex- 
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ternal stationary monochromatic field. The dipoles are 
spheres of radius a — 0.1/im and real dielectric constant 
e e = 16. For this model, the bare electrostatic polariz- 
ability is a = An^^a 3 ~ 1.11 • 10~ 20 m 3 . The wave- 
lengths of the probe field range from 300nm to 500nm 
which correspond to frequencies between 3 • 10 5 m -1 c and 
5 • 10 5 m~ 1 c, much lower than any resonant frequency by 
assumption. To ensure that both the dipole approxima- 
tion and the low frequency condition hold, the radius 
of exclusion volume is set to £ = OAfim. In order to 
avoid overdensity terms in the correlation function and, 
at the same time, allow for multiple-scattering, we set 
(\> = 0.4. Note that because £ = 4a, the material fill- 
ing fraction is 8 times less, 0.05. This corresponds to a 
numerical density p = 1.49 • 10 19 m" 3 . We perform the 
computation of the refractive index and the extinction 
coefficient in three approximations. In the first one, a 
non-self-consistent MG approximation is used. In it, sin- 
gle particle polarizabilities are only affected by in-free- 
space radiative corrections, a = ocq\1 — g^A^ao] -1 , and 
the effective susceptibility adjusts to that of the s.s. MG 
model. In the second one, the self-consistent overlap MG 
approximation is employed and the 7-factors computed 
in Section IIIII are used to renormalize both the single- 
particle polarizability and the dielectric constant. In the 
third one, the self-consistent overlap recurrent approxi- 
mation of Section [V] is used. The results are plotted in 
the graphs of Figj5] 

We have also computed the total power emitted 
by a random dipole when stimulated by a stationary 
monochromatic field in the same range of frequencies. 
Eq. (l65l) is applied for W Tot \ For the computation of 
the coherent emission, W Coh \ we use Eqs. (|45l64[) . The 
results are plotted in the graphs of FigJU The val- 
ues are normalized to that of the in-free-space emission, 
W° = |e 3{a}. 



VII. DISCUSSION 

A. Corrections at high filling fractions 

When the effective filling fraction of scatterers is rela- 
tively high, the first-neighbor overdensity term in h(r) 
becomes relevant for the computation of x^ ■ That 
is the case when the relation £/p -1 / 3 approaches one. 
The functional form of the overdensity term depends on 
the precise profile of the interaction potential between 
dipoles. For the sake of simplicity we will restrict to the 
delta function of Eq.(p]l. Ah(r) = C £6 {1} (r - £). The 
inclusion of such a term in the computation of x^ (Q) 
yields the additional contributions, 

A X i 2) (Q) * -C{ P af[\e + ^Q 2 + 0(C 3 X 2 Q 2 ) + ■■], 

(70) 




FIG. 5. Refractive index, n, and extinction coefficient, 
K, for the numerical example of Section IVI CI The long- 
dashed line represents the non-self-consistent MG approxi- 
mation (nScMG). The short-dashed line corresponds to the 
self-consistent overlap MG approximation (ScOvMG). The 
solid line represents the self-consistent overlap recurrent ap- 
proximation (ScOvR) of Section [V] The curves of n appear 
separated in the lower plots. 



( 71 ) 

The relevance of these terms in our previous formulae 
is as follows. Regarding far-field calculations and co- 
herent emission terms, only long wavelength modes with 
Q <C 1 enter Eqs. (|26l58[) . Those equations contain the 
order C(C°) terms which amount to the transverse coher- 
ent modes. The inclusion of Eqs. (l7QI71[) is irrelevant for 
them. Therefore, the degree of accuracy of Eqs. (|26l58p 
remains unaltered regardless of the degree of order in the 
medium. The same argument applies to Eqs. (l45l64j) . On 
the contrary, Eqs. (j27l33|) and Eqs. (|59H62j) contain short 
wavelength modes with Q > 1. Therefore, in order to 
keep their degree of accuracy for high filling fractions, 
it is necessary the inclusion of additional terms. Those 
are, the ones of order <D{(°) of Eq.® in Xxmg™ > Xx 5 
and those of orders 0(£°), C(C 2 ) and C(C 3 ) of Eq.(|7T]l 

(2) -(2) 
m X||MG°"'^|| ' 

Regarding experimental observations, in order to pre- 
vent high-order effects in the artificial solids, it would 
be necessary to deform the surface of the silica shells 
to avoid first-neighbor overdensities around each coated 
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FIG. 6. Plots of the total and coherent stimulated emission 
in the self-consistent overlap MG approximation (ScOvMG, 
short-dashed line) and in the self-consistent overlap recurrent 
approximation (ScOvR, solid line). Quantities are normalized 
to the in-free-space emission, W°. 



gold particle. 



B. Comparison with previous approaches 

It has been mentioned already that numerical simula- 
tions do not take proper account of radiative corrections 
in the renormalization of the single-particle polarizabil- 
ity. Hence, the dipole emission cannot be evaluated by 
direct application of the optical theorem in complex me- 
dia, Eq. (|f)5i) , unless knowledge of Xeff be implicitly as- 
sumed in some approximation -eg. Ref. HH Nonethe- 
less, we have explained how both the total and coherent 
intensities can be evaluated in simulations by exciting 
one of the dipoles of a sample with an external probe 
field as in Ref. [28| and |46|. Mallet et aZ.— have measured 
instead the coherent and incoherent intensity scattered 
by a finite sample. They have compared the numerical 
data with the corresponding formulae for a Mie sphere 
with an effective dielectric constant adjusted to the MG 
formula. Our reasoning in Section IVII Al would explain 
why the authors found good agreement for the coherent 
cross-section, even for high filling factors. On the con- 
trary, because the ordinary MG formula does not include 
the C-dependent terms which enter the incoherent inten- 
sity (or the macroscopic absorption in the terminology of 
Ref. |23[ ) , it is not surprising the large discrepancy found 
also there in the incoherent cross-section for high effec- 
tive dielectric constant. Direct comparison of our predic- 
tions with the numerical data of Ref. [23| is not possible 
as the authors deal with scattering by a Mie particle. It 
is however possible to verify that, in the single-scattering 
approximation and in the small particle limit, our de- 
composition between coherent and incoherent emission 
agrees with the calculation of Froufe et a/.—. 

Cichocki and Felderhof developed in Ref. IH a stochas- 
tic approach to the dielectric constant which is formally 
very similar to ours. Our self-polarization 7-factors of 



Eqs. (|7l6[) are analogous -but not physically equivalent - 
to their reaction field operator. Our two-scatterer suscep- 
tibility in the quasicrystalline approximation, x^ 2 \ 1S an 
approximation to their short-range connector. As a re- 
sult, they obtained a formula for the dielectric constant 
which is of the same functional form as that in Eq. fTTj]) 
once the limit q — > is taken. A similar strategy was 
followed by Leegwater and Mukamel in Ref. [26|. That is, 
they first dressed up the single particle polarizability in- 
cluding near-field interactions and then applied the MG 
formula over the renormalized polarizability. Cichocki 
and Felderhof applied their formalism to the particu- 
lar case of the Drude-Lorentz model (equivalently, the 
Maxwell- Garnett model) in a non-polar flui d 12 i 13 ' 32 . In 
the two latter references they concentrated on the reso- 
nant regime. For the computation of the reaction field 
operator they developed in Ref. [H a self-consistent ring 
approximation which is similar to our self-consistent over- 
lap recurrent approximation. However, our approxima- 
tion is performed off resonance and for the sake of sim- 
plicity we restrict to 0[a 2 /(47r£ 3 ) 2 ] terms in x^ rec '. Our 
approach includes transverse modes while the authors of 
Ref. [H, HH, and |H restricted to electrostatic induction 
only. More importantly, because transverse modes are 
included neither in Ref. [H, [ID, and HH nor in Ref. HH, 
the radiative modes of Eqs. (|26l58l) are missing there and 
the computation of the radiative line-width cannot be 
complete. As a matter of fact, neither a proper renor- 
malization of the single-particle polarizability that satis- 
fies the optical theorem nor the formula for W can 
be obtained following those approaches. Also, while the 
zero- modes which enter our Xeff are properly accounted 
for in Ref. [l2|, [l3|, and [HI since they are mainly affected 
by electrostatic interactions, only the 0((~ 3 ) terms in 
our 7|| can be obtained out of those works. 

There are two fundamental differences between our ap- 
proach and that of Cichocki and Felderhof. The first one 
is technical. Because we have obtained in Ref. [l] exact 
formulae for Q(r, r) as a function of x{q) -Eqs. (|3l4p . our 
strategy consists of computing first x(q) i n some sensi- 
tive approximation. While its zero-mode is by definition 
Xeff, application of the exact formula for G(f,r) yields 
the radiative corrections which enter a . On the contrary, 
Cichocki and Felderhof developed in Ref. |l_2j an approx- 
imation to capture the most relevant recurrent scatter- 
ing diagrams for the computation of a. Those diagrams 
account for the multiple reflections that actual probe 
photons experience within clusters of particles, restricted 
to electrostatic induction. In their approximation, con- 
secutive scattering processes are never uncorrelated and 
hence Dyson's propagator is never present. Because it is 
Dyson's propagator the relevant one for radiation, that 
approximation could not capture all the radiative cor- 
rections, but for those in free space at the most, even if 
radiative induction were implemented in their approach. 
In Ref. 13 the authors developed another approximation 
for the computation of their connector whose zero-mode 
enters Xeff- To this respect, ours is a lower order ap- 
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proximation. 

The second difference is conceptual and more subtle. 
The self-consistency requirement that the ring approxi- 
mations of Ref. [HI and [l3l satisfy are in perfect agreement 
with the application of the Couple-Dipole-Method. In 
fact, self-consistency is inherent to the CDM equations. 
However, the self-consistency that we demand on radia- 
tive corrections apply to virtual photons instead. Hence, 
our comment at the beginning of Section IVI1 That is, 
it is recognized that the application of the CDM equa- 
tions needs of some a priori prescription for the value 
of the radiatively-corrected single-particle polarizability 
in order to obtain physically acceptable values for the 
scattering cross-section -see eg. Ref. HH. Our approach 
reveals that self-consistency in radiative corrections re- 
quires to solve for a at the same time that the CDM 
system of equations is solved^. 

Nevertheless, the elegant analytical approach devel- 
oped in Ref. [l^ andfTijl will pave the way for the computa- 
tion of the radiative corrections at resonance in a future 
publication. 

C. Relationship between the emission spectrum and the 
local field factors 

In Section TVI Bl it has been emphasized the nearly uni- 
versal character of Eqs. (l45l66|) and Eqs. (|64l69l) in the low 
frequency limit of the virtual cavity scenario. Related 
with this fact, there exists a strong motivation to vali- 
date Eqs. (|45l66|) and Eqs. (|64l69|) . That is to show that 
the formulae for coherent emission in the stricto sensu 
virtual cavity scenario contain only one local field fac- 
tor, Sft{ x± ^ — -}, and not two as commonly claimed in 
the literature -eg. Ref. [H [M HHH In the follow- 
ing we will restrict ourselves to the s.s. MG model. 

For this model the relations = SR(eM 3 G}+2 , 

ijftj *mg 1 _ ^'mcH 1 are exac t i n three and two dimen- 

sions respectively and so are Eqs. (|45l66l) . Therefore, we 
have proved explicitly that only one Lorentz-Lorenz local 
field factor appears in the formula for coherent emission 
in an MG dielectric. 

In the literature, the usual expression in a three- 
dimensional dielectric for which an effective dielectric 
constant can be defined reads, 

W^cxp^) 2 ^, (72) 

where absorbtion is ignored. The terminology varies 
from one author to another, referring to the above for- 
mula as virtual cavity 5 -, Lorentz-LorenaS or Claussius- 
Mossotti 5 ^ emission^- Certainly none of those pioneer 



Most of those papers deal with the spontaneous decay rate in- 
stead. The distinction is irrelevant for the present discussion. 



physicist worked out Eq. (l72l but just the relation be- 
tween macroscopical classical fields in effective media* 2 - - — . 

It is clear that, if the above equation refers to the 
s.s. virtual cavity scenario and an effective medium ex- 
ists, the medium must be an MG dielectric to which 
all the formalism developed throughout this work ap- 
plies. In the papers cited above it is not clarified the 
coherent/incoherent nature of the radiation in Eq. (l72l) . 
However, the presence of the refractive index factor, 
n = sje^j], denotes the integration of J 2^{Gj_(q)} 
which would yield only coherent power. Thus, the ex- 
pression in Eq. ([72]) should stand for coherent emission. 
In any case, because the transverse incoherent emission in 
our calculation contains only terms of order 0(£ _1 ), the 
addition of transverse coherent and transverse incoher- 
ent emission cannot equal Eq. ([72l) . Also contrary to the 
general assumption, effective longitudinal modes carry 
incoherent radiation even in the absence of absorption 27 . 
That is proportional to the C(C°) terms of 7* /G ""' . The 
reason why some experiments on spontaneous emission 
seem to agree with Eq. (|T2"|) 60 is that those experiments 
measure total far-field radiation, and not only either co- 
herent or transverse. It turns out that for n close to 
one Eq. (|?2")) approximates the value of the total far-field 
radiation. That is the sum of Eq. (l4"5|) plus longitudinal 
incoherent radiation. This is illustrated in FigJT] where 
the data refer to the numerical example of Section IVI CI 
Therefore, the rule out of Eq. d72l in favor of our result 
could be verified by direct measurement of the coher- 
ent power. Note also that corrections due to recurrent 
scattering are of order O(e^j) and higher. Nevertheless, 
Eq. (IT2")) fails already at leading order. 

We stress that our conclusion about the presence of 
a unique local field factor in Eq. (|45|) is not model- 
dependent. The equivalence jft{ Xx( ^°' ] } = ^ e ™^+ 2 
is of course only valid for the s.s. MG model since it is 
implicit that both the condition fc£ <C 1 holds and the 
exclusion volume function is the only correlation func- 
tion. Nonetheless, Eq. (|4"2"j) is exact in the s.s. virtual 
cavity scenario. That is, if the emitter can be treated as 
a single dipole which is in all equivalent to the rest of 
the constituents of the medium. As explained in Ref. 1. 
and [60l . it is only for stimulated emission that the na- 
ture of the emitter remains equal to that of the rest of 
dipoles during the emission process. Hence, our proposal 
to measure W ( F oh - instead of T CJoh - . 



VIII. CONCLUSIONS 

We have computed the radiative corrections to the off- 
resonance dielectric constant of an MG dielectric. To 
that aim, we have used the formalism developed in Ref. QJ 
As a byproduct, we have obtained formulae for the spec- 
trum of the coherent stimulated emission in both two 
and three dimensions. Two variants of the quasicrys- 
talline approximation have been used. The first one is 
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Appendix A: Coherent emission of TE modes in two 
dimensions 

In this section we derive Eq. (|66[) for a two-dimensional 
dielectric made of cylindrical inclusions homogeneously 
distributed as in FigfJJ Cylinders of radius a -C A and 
permitivity e e have a bare classical and isotropic polar- 
izability a° = 2 na 2 j^j for p-polarized field o 51 ' 61 . The 
Maxwell-Garnett formula reads, 



FIG. 7. Plots of the total and coherent stimulated emission 
for the numerical example of Section TVI CI The solid line cor- 
responds to the total emission in the overlap self-consistent 
MG approximation. The long-dashed line represents the 
transverse emission according to Eq. (|72[l . The short-dashed 
line corresponds to the coherent emission in the overlap self- 
consistent MG approximation. Quantities are normalized to 
the in- free-space emission, W°. 



referred to as the self-consistent overlap MG approxima- 
tion. In it, the only relevant correlation between dipoles 
is approximated by a spherical exclusion volume. The ef- 
fective susceptibility is that of the ordinary MG formula, 
Eq. (|16p . The 7- factors which enter the renormalization 
of a are those of Eqs . • The formulae for coherent 

radiation are those of Eq. (|45|) and Eq. ([66|) in three and 
two dimensions respectively. The second approximation 
is referred to as the self-consistent overlap recurrent ap- 
proximation. In it, in addition to the spherical exclusion 
volume, self-correlation is also considered and recurrent 
scattering diagrams are included up to 0(a 2 /(47r£ 3 ) 2 ) in 

Xst c ' '. The effective susceptibility is that of Eq.(f56"|). 
The 7-factors which enter the renormalization of a are 
those of Eqs. (|58H62l) . The formulae for coherent radia- 
tion are those of Eq. tjM]) and Eq. (|6T)]) in three and two 
dimensions respectively. In Section IVT1 we have proposed 
several numerical and laboratory experiments to test our 
results. Most of them are based on the measurement 
of the coherent stimulated emission since the functional 
form of its formulae is nearly independent of the micro- 
scopical details of the medium. To this respect, we have 
found that only one local field factor enters those formu- 
lae in contradiction with the common assumption that 
two factors do. 
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so that the local field factor is 



(Al) 



(A2) 



The dipole moment induced by an external field E^r) 
parallel to the plane on the emitter rod at r is p(r) = 
eoctplEg (r), I being the length of the rod. The 
Dyson propagator of transverse TE modes in the long- 
wavelength limit is, 



Gi D (q,u) 



-MG 



k 2 



(A3) 



Correspondingly, the spectrum of coherent bulk modes 
reads, 



3{ 7 i D } buife = 3{ 



d 2 g 

(27T) 2 



1 + 7r arctan 



In function of Eqs. (|A2IA4| , the coherent power of TE 
modes is, 



(A5) 



4c^ aa p 



Substitution of Eqs. (|A2IA4p in the expression above 
yields Eq.([66]). 
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